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Abstract 

We study the effects of random bonds on spin chains that have an excitation 
gap in the absence of randomness. The dimerized spin- 1/2 chain is our princi- 
pal example. Using an asymptotically exact real space decimation renormal- 
ization group procedure, we find that dimerization is a relevant perturbation 
at the random singlet fixed point. For weak dimerization, the dimerized chain 
is in a Griffiths phase with short range spin-spin correlations and a divergent 
susceptibility. The string topological order, however, is not destroyed by bond 
randomness and dimerization is stabilized by the confinement of topological 
defects. We conjecture that random integer spin chains in the Haldane phase 

exhibit similar thermodynamic and topological properties. 
Pacs: 75.10.J, 75.30.H, 75.50.E 
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Extensive theoretical work on random quantum magnetic systems has been carried out 
since the late 1970's JI]-|§]. Systems that behave critically in the absence of randomness are 
unstable against weak randomness and flow to the random singlet (RS) phase JI|||. In the 
RS phase, spins far apart in space form weakly bound singlet pairs in a more or less random 
manner. This phase is also referred to as the valence bond glass (VBG) phase 0. The 
low temperature thermodynamic properties of these systems are dominated by the weakly 
bound pairs and are universal For instance, the susceptibility of the undimerized 

random antiferromagnetic Heisenberg or XX Z spin-1/2 chain diverges as [Tlog 2 T]- 1 at low 
T, independent of the details of the randomness ||. Universal power law behavior has also 
been found in disorder averaged spin-spin correlation functions ||. Experiments, however, 
seem to find power law divergent susceptibilities with nonuniversal exponents 0. It is of 
interest to study if there exist relevant perturbations at the RS fixed point that drive the 
system towards a state exhibiting the nonuniversal behavior found experimentally 

A related issue is the effect of randomness on spin chains that have an excitation gap in 
the absence of randomness || . The most prominent examples of such chains are integer spin 
chains in the Haldane phase |5J. Other examples include dimerized spin-1/2 chains [|IU| and 



spin chains with spontaneous dimerization JTT| , [n^] . All of these systems have topolologically 
ordered |13|] ground states. One might think that strong enough randomness will inevitably 



destroy the topological order of the ground state. However Haldane has suggested |L4| that 



there exists a class of random perturbations for which the the topological order in the ground 
state of integer spin chains is stable regardless of the strength of the perturbations. We will 
show that the analog of this prediction for random bond dimerized spin-one-half chains is, 
in fact, correct. 

An explicit example that provides strong support to Haldane's conjecture (sic!) can be 



found in the random version of the AKLT model: [ 15 ] 

1 



Si • Si + i + —(Sj • Si + i) 2 



(1) 



where Jj > 0. The exact ground state of this random model is identical to that of the pure 



model, [|15[ i.e., a valence bond solid. Its excitation spectrum and thermodynamic properties 



will certainly depend on the distribution of J i; yet the perfectly topologically ordered ground 
state is completely unaffected by randomness. 

In the absence of disorder, spin-one chains in the Haldane phase and dimerized spin-one- 
half chains exhibit similar physical properties |16|| . They both have a nondegenerate ground 
state with an excitation gap, and more importantly, they both have string-topological order 
16fl . Hida has shown that they can be continuously connected to each other without closing 



the gap or removing the topological order; i.e., they are in the same phase. [|16| It is natural 
to expect that they also behave similarly in the presence of randomness. 

In this paper we study the random bond dimerized spin-1/2 chain in detail. Using 
the asymptotically exact real space decimation renormalization group introduced by Ma, 
Dasgupta and Hu and extended by Fisher ||, we find that enforced dimerization is a 
relevant operator at the RS fixed point that drives the system to a random dimer (RD) phase. 
The low temperature thermodynamic properties of the RD phase are nonuniversal. For 
weak dimerization, the spectrum of the RD phase is gapless and the susceptibility diverges 
as x ~ T~ 1+Q with a > and dependent on the bond distribution, (similar to the behavior 
found experimentally and in qualitative agreement with the RS thermodynamics), but the 
averaged spin-spin correlation function remains short ranged. Thus for weak dimerization 
the RD phase is an example of a Griffiths phase. More importantly, we find that the string- 
topological order is not destroyed by random bonds. We conjecture that these results also 
apply to random bond integer spin chains in the Haldane gapped phase. Comparison will 
also be made with spontaneously dimerized spin-chains which behave very differently upon 
introducing disorder. 

Consider the model Hamiltonian 

H = Ji [Sf Sf +1 + SfSf +1 + AS*Si +1 , (2) 

i 

where are spin-one-half operators, Ji are (random) positive coupling constants, and 
< A < 1. Here we will concentrate on the cases A = (XX chain [|17]) and A = 1 



(Heisenberg chain), since it has been shown that, for the case of random bonds, the Ising 
coupling is irrelevant when A < 1 and that the system flows to the random XX chain in 
the low energy limit ||. We assume that the distribution of the couplings, Jj, depends 
on whether i is even or odd. We write the distribution functions for even and odd J's 
as P e (J, J ) and P Q (J, J ) respectively. Here J is the cutoff in the distribution function 
corresponding to the strongest bond in the system. As Fisher |J has shown, in the absence 
of dimerization, i.e., when P e (J, J ) = Po(J,Jo), the low-energy, long-distance behavior of 
Eq. (0) is universal, and the XX and Heisenberg chains behave in essentially the same way. 

Following Fisher ||, we introduce a decimation renormalization group procedure, in 
which we pick the bond in the system with the largest J (i.e., the strongest bond), say J 2 
between spins 2 and 3. Since this is such a strong bond, spins 2 and 3 are likely to form a 
singlet pair and become unimportant at low energies (on scales much smaller than J 2 ). The 
major physical effect of the existence of spins 2 and 3 is to generate an induced coupling 
between their neighboring spins 1 and 4. For the XX chain: -f/i-4 = Ju(S*S% + SfSf) 
where J u = J1J3/J2 + 0(1/ Jf) and for the Heisenberg chain: = J14S1 • S 4 where 

J14 = J 1 J 3/ '(2 J 2) + 0(1/ Jf). The effect of this decimation procedure is to get rid of the 
strongest bond (and also its two neighbors) in the system, generate a weaker bond between 
the spins neighboring the decimated ones, and lower the overall energy scale. This procedure 
becomes asymptotically exact in the low energy limit ||. The new energy cutoff is then 
lowered to Q = max{J}. Following Fisher and anticipating that the bond distribution will 
become broad on logarithmic scales at low energy || , we transform to logarithmic variables 
and define T = — Iog(f2/J ) and ( = log(f2/J)/r, so that both T and ( are positive and a 
larger T and a larger ( correspond to a lower energy scale and a weaker bond respectively. 
The recursion relations now become (keeping the leading term only) 

Ci_ 4 = Ci + Cs - C2 + « = Cl + Cs + « (3) 

where we used the fact that ( 2 — since J 2 = is the strongest bond in the system, k = 
for the spin-one-half XX chain and log(2) for the spin-one-half Heisenberg chain. The flow 
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equations for the bond strength distribution functions p e ((,Q) and p ((,Q) in terms of ( 
are then 

<9p e ,o(C,r) <9p e , 

— — = -q^- + [p e ,o(o, r) - p 0>e (o, r)]p £i0 

+ Po , e (o, r) J J d&dCzp^id, r)p e , (c 2 , r)5(c - d - C2 - «)■ (4) 

When k = (i.e., the XX chain case), these flow equations are identical to those encountered 
in the transverse field Ising model if we identify the even bonds as the bonds between Ising 
spins and odd bonds as the transverse fields |HJI§fl- 111 order to find fixed point solutions of the 
renormalization group (RG) flow, it is necessary to rescale variables. Following Fisher [|5]||, 
we introduce the rescaled variable r\ = (/T and the new distribution function Q ej0 (r],T) = 
rP e o ( J, Q). The flow equations for Q are 

= Q et0 + (1 + rf}^ + [Q e , o (0, r) - Q o>e (0, r)]Q e , 

+Q Oi e(0, r) J J d?7iCf%Qe, (^l)Qe, (%)<% ~ Vl ~ V2 ~ f) ■ (5) 



As Fisher has shown ||[5|], the flow equations (^) have only one generic fixed point [HJ 



Qe = Qo = Q*(v) = e~ v ®(v)- (6) 

This fixed point distribution corresponds to the random spin-1/2 chain without dimerization, 
a model studied extensively before Going back to the original variable (, we find the 
fixed point distribution corresponds to 

P(0 = fe~^, (7) 

i.e., the width of the distribution on the logarithmic scale grows linearly with the log of 
the energy scale T. For small deviation away from the fixed point Q e = Q* + q e and Q Q = 
Q* + q Q , there is only one relevant eigenperturbation behaving as q et0 (r],T) = q et0 (r])T x 
with eigenvalue A = 1, and the eigenvector is q e = (rj — l)e~~ v and q a = —{t] — l)e~ v . The 
relevant perturbation, like the fixed point distribution, is independent of k; hence the XX 
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and Heisenberg chains will behave similarly. However, non-zero k does affect the irrelevant 
perturbations. 

The relevant perturbation corresponds to the difference in the distributions for even and 
odd bonds. Therefore we find that dimerization is a relevant perturbation near the RS fixed 
point, with eigenvalue +1. 

For weak dimerization, the system barely knows that there is a small difference between 
even and odd bonds in the early stages of the RG flow. Both distributions initially flow 
toward the RS fixed point solution with a small relevant perturbation reflecting the extent 
of dimerization: 

Qo(T) =Q* + 6T(r l -l)e- r >, 

Q e (T)=Q*-5T(r,-l)e- r ' (8) 

where 5 characterizes the strength of the dimerization (distance from criticality). In general, 
5 depends in a complicated way on the shape of the original distributions and at what energy 
scale it is defined. As the flow away from the RS point continues, the even (odd) bonds 
get much weaker than the odd (even) bonds if originally the even (odd) bonds were only 
slightly weaker than the odd (even) bonds. The relevant perturbation grows linearly with 
T and becomes of order 0(1) as T = T ~ 1/|5|. The flow equation for the density of spins 
that have not yet formed singlets at energy scale T is 

^ = -2«(o,r) (9) 

Using the RS fixed point distribution Eq. (H) in Eq. (|9|) we find that n ~ 1/T 2 [|,|J so that 
when the relevant perturbation becomes large, the density of active spins is n ~ 5 2 . The 
corresponding length scale L, which is the typical distance between the remaining spins, is 
L ~ Tq ~ 1/5 2 . At this stage the existence of dimerization becomes dominant and under 
RG most of the bonds decimated are odd bonds. 

The fact that a small difference in the bond distributions grows as one lowers the energy 
is physically easy to see. Assume the odd bonds are slightly stronger than the even bonds 



in general. Then in the decimation procedure, it is slightly more likely that an odd bond 
gets decimated. When that happens, typically two intermediate strength neighboring even 
bonds also disappear, and a much weaker even bond gets generated. Hence, the width of the 
even bond distribution grows faster than the odd bond distribution, and its over all strength 
also decreases faster. Thus, in the low energy limit, the system can be viewed as a trivially 
soluble collection of uncoupled spin pairs (isolated odd bonds). We refer to this phase as 
the Random Dimer (RD) phase. 

After renormalization the distribution of odd bonds takes the form 

Po(0 ~ ^e-^e(c). (io) 

1 o 

In terms of the original variables the odd bond distribution is 

Po(j) = 7f(7fr 1+Q ©(i - <f ), (ii) 

where Qo = Joexp(— r ) and a oc 5. 

This effective independent pair Hamiltonian with a power law bond distribution is iden- 
tical to that introduced by Clark and Tippie |2D| to explain the low temperature thermody- 



namics of the random spin chains. Here we have derived it using RG from a realistic model. 
The leading temperature dependence of thermodynamic properties can be determined by 
assuming that all spins connected by bonds with energy greater than the temperature have 
paired up into singlets and all spins connected by bonds with energy less than the tem- 
perature are essentially free. This is a good approximation for broad bond distributions. 
In this way, the specific heat and susceptibility in the low temperature limit can be easily 
calculated. As the temperature goes to zero, the the spin susceptibility (in any direction, 
with possible direction dependent prefactors) diverges like \ ~ T"" 1 , and the specific heat 
goes to zero like C v ~ T a . The averaged spin-spin correlation function is short ranged, 
with the correlation length (distance between spins) £ ~ \8\~ u ~ |^|~ 2 - The existence of a 
divergent magnetic susceptibility away from the critical point is characteristic of a Griffiths 



phase. The divergent susceptibility arises from magnetically active gapless excitations. [21 



For the RS phase discussed by Fisher |5j , the averaged spin-spin correlation function decays 
as 1/R 2 at long distance so the system is critical and one expects a divergent susceptibility 
The Griffiths phase in the random dimerized spin-one-half chain is exactly analogous to the 
Griffiths phase that appears in the random transverse field Ising chain ||. If the initial 
dimerization is large, the flow begins far from the RS phase and the bond distribution of 
the stronger bonds does not flow to a power law and the gap does not close up. When this 
is the case, thermodynamic properties will depend strongly on the initial distribution and 
the susceptibility will remain finite. 

The dimer phase has a novel kind of topological order that measures the dimerization of 
the chain. The "string-topological correlation function" is 



T 



ij 



SI exp 


^ E si 


s; 




i<k<j 





(12) 



where \^q) is the ground state. Ty is similar to the topological correlation function for the 



spin-one chain |L3[ and it maps onto the spin-spin correlation function of the transverse field 

-1/4 if i is a left spin of 



Ising chain. |22| For a completely dimerized ground state - 
a dimer and j is the right spin of a (possibly different) dimer. This is because every spin 
between i and j in the completely dimerized model is paired up with another spin between 
i and j. Tij = otherwise. Therefore this special topological correlation function is long- 
ranged although there is only short range spin-spin correlation, a situation similar to the 
special kind of off-diagonal long range order (ODLRO) in the fractional quantum Hall effect 
(FQHE) p3[] . I n the pure case this topological order is closely related to the existence of a 
gap |T3|] , again just like in the FQHE. 

For a general spin-one-half chain with randomness, we introduce a topological order 
parameter 



T = lim (T 2i2 i-i 



To,; 



2i+l,2j), 



(13) 



where the overbar stands for average over randomness. In the absence of dimerization, T 
vanishes. For a random system, T measures the probability that the two end spins survive 



S 



decimation until the dimerization becomes large, and the low energy physics becomes that 
of the completely dimerized chain. This probability is just the square of the density of spins 
at the dimerization crossover scale. Therefore, for small 5, T scales like 

T ~ -\6\^sga(6), (14) 

with (3 = u = 2, and sgn(<5) is + if even bonds are stronger and — if odd bonds are stronger. 
In the absence of randomness f3 is found to be 1/12 . 

The dimerization described in this paper is premeditated because its magnitude and sign 
is determined by the details of the bond distributions. If the even bonds are stronger than 
the odd bonds then it costs less energy to form singlet pairs over even bonds than over 
odd bonds. In the pure case, premeditated dimerization leads to a unique ground state. 
In contrast, spontaneously dimerized spin chains can dimerize over even or odd bonds with 
equal energy cost. In the pure case, spontaneous dimerization leads to degenerate ground 
states. The fact that there is a right way and a wrong way to dimerize a premeditated 
chain and no wrong way to dimerize a spontaneous chain leads to profound differences in 
the response of these systems to random perturbations. 

For the case of premeditated spin chains, the random perturbation and quantum and 
statistical fluctuations will generate regions where the chain has dimerized in the wrong way. 
Between regions of even and odd dimerization there has to be a topological soliton which 
is essentially a spin unpaired to its neighbors. In the presence of dimerization, the energy 
cost of the wrong region is proportional to its length (at least at long enough length scales); 
hence the topological solitons are confined by a linear confining potential and the topological 
order persists |24[]. The confinement length in the weak dimerization limit is essentially the 



length scale at which dimerization becomes significant under RG, which is also the spin- 
spin correlation length. We hence find that although the gap vanishes in the presence of 
strong randomness, the dimer phase is stable and the topological order persists due to the 
confinement of topological defects. The situation for the FQHE is similar: incompressibility 
is lost upon introducing randomness, while ODLRO survives [P5fl . Another example of this 
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phenomena occurs in dirty superconductors. Moderate disorder can remove the energy gap 
without destroying the condensate. 

In contrast, topological solitons are unconfined in spontaneously dimerized systems be- 
cause there are no wrong regions. Therefore, unlike the case of premeditated order, weak 
randomness will destroy spontaneously generated dimerization, despite of the existence of a 
gap in the absence of randomness. 

We have shown that the dimerized spin-one-half chain is stable against disorder. This 
can be seen as due to the confinement of topological defects. When the dimerization is 
weak or randomness is strong the system is in a Griffiths phase in which the spin-spin 
correlation function is short ranged yet the susceptibility diverges. The susceptibility and 
specific heat follow simple power laws with nonuniversal exponents at low temperatures. The 
string-topological order is not be destroyed by disorder. We conjecture that these results 
apply to other systems with premeditated (not spontaneously generated) topological order 
such as the spin-one chain in the Haldane phase. We also conjecture that spontaneously 
generated topological order is unstable against disorder. A detailed analysis using a real 
space decimation procedure that is appropriate for the spin-1 and higher spin chains will be 
presented elsewhere Experimentally there has been some work on effects of hole doping 
in the spin-1 chains f27fl- The effects of hole doping is probably very different from random 
bonds in both spin-1 chains and dimerized spin-1/2 chains because it introduces unconfined 
topological defects into the system. Further investigation is underway . 

We are extremely grateful to Daniel Fisher and Duncan Haldane for numerous stim- 
ulating discussions and important suggestions. This work was supported by NSF grants 
DMR-9224077 and DMR-9416906. 
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